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Abstract 

We consider the hydrodynamic regime of gauge theories with general triangle anomalies, where 
the participating currents may be global or gauged, abelian or non-abelian. We generalize the 
argument of arXiv: 0906. 5044, and construct at the viscous order the stress-energy tensor, the 
charge currents and the entropy current. 
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I. INTRODUCTION AND SUMMARY 



In this letter we will consider the hydrodynamic regime of microscopic quantum gauge 
theories with triangle anomalies. The field theory currents can be global or gauged, abelian 
or non-abelian. We will label this full set of currents with the indices (a, /?,...), and will 
limit ourselves to symmetries which are not spontaneously broken. 

Some of the field theory charges will be accessible at the energy scale defined by the 
temperature. Such charges and their chemical potentials may take part in the hydrodynamic 
description, as part of the local equilibrium parameters. There are possible obstructions to 
this, e.g. a gauged charge may be screened. Also, to be part of the thermodynamics a charge 
must be conserved with sufficient accuracy within each local equilibration region; thus, a 
global U(l) charge will not be part of the fluid description if its conservation is too disrupted 
by anomalies with gauged currents. For the bulk of the discussion, we will assume that the 
conservation of hydrodynamic charges is not violated at all. However, during the derivation 
we will violate the conservation of charges using fictitious external fields. The more general 
case of slight violation will be discussed in section HVl 

Gauge fields may also enter the fluid description by assuming nonzero average values over 
macroscopic distances. This possibility upgrades the discussion from pure hydrodynamics to 
electro/magneto-hydrodynamics. If a gauge field enters the description, so will the current 
that it induces; the associated charge density may also appear as an independent parameter, 
or it may be screened. 

For non-abelian gauge fields and their associated charges, a complication arises. If a 
charge density is to be part of the thermal parameters, it must be additive and conserved 
(in the ordinary, not in the gauge-covariant sense) within each local equilibration region. 
In other words, within each equilibration region the charge must be approximately abelian. 
This imposes a condition on the fluctuating microscopic component of the associated gauge 
field. The same condition must hold for the gauge field itself to enter a fluid description: 
without a consistent gauge frame within each equilibration region, the gauge field cannot 
have a smooth, well-defined average value. The technical consequences of this condition are 
that gauge-covariant derivatives commute with thermal averaging, and that the averaged 
field strength can be derived from the averaged potential A a ^. 

Our aim is to find the allowed leading viscous order constitutive relations of the field 
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theory hydrodynamics. In order to keep track of orders of magnitude, we introduce a formal 
small parameter e. The "order" of a quantity will refer to the power of e involved. We will 
take all the gradients and field strengths to be of order ~ e. The background metric g^ v 
may be curved on the length scale of the hydrodynamic gradients, i.e. its Riemann tensor 
is of order 0(s 2 ). 

The zeroth-order constitutive relations for the stress-energy density, the charge currents 
and the entropy current have the standard ideal-fluid form: 

J { a^ = V=gn a u» , (1) 

where is the energy velocity, = — 1, e is the energy density, p is the pressure, n a is 
the charge density and s is the entropy density. P" — 8 V + u^u v is the projector orthogonal 
to u^. The indices (a, b, ... ) enumerate the currents which participate in the hydrodynamics 
(whether through a charge density or a gauge field). For screened charges, we set n a = 0. T£ 
is the "intrinsic" stress-energy density, i.e. it includes the contribution from the microscopic 
component of the gauge fields, but not the contribution from the macroscopic gauge fields. 

Our result for the first-order constitutive relations is a direct generalization of [1]. We 
will find (in the Landau frame): 

4^ = v^g° a b (K - TP^D U ^) + ^ + eifs 6 " , ( 2 ) 

where T = (de/ds) n is the temperature, [i a = (de/dn a ) s are the chemical potentials (with 
Ha = for screened charges) , r\ > is the shear viscosity, ( > is the bulk viscosity, and a a b 
is the conductivity matrix with a^ ao ) positive semi-definite. is the covariant derivative 
with respect to the metric and the gauge fields. The derivative of \i a jT is taken in the 
gauge-covariant sense, i.e. D v \i a = d u fi a + f a b c A b u fi c , where / a & c are the Lie algebra structure 
constants. The shear tensor ir^ u and the vorticity density w M are defined by: 

V, = P£P°D {p u a) - \p^D p u p , 

! 3 (3) 

^ = _ e w UvdpU(j ; 
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where e^ upa is the metric- independent Levi-Civita density with components ±1. The electric 
field E£ and magnetic field B£ are defined by: 

E^ = Fru v ; B^ = ^V^, (4) 

where F* is the macroscopic gauge field strength. For global charges, we set = 0. 

The ordinary viscous coefficients (, rj and a a b can be arbitrary functions of state. In 
contrast, the vorticity/magnetic coefficients £ a , £ and £i B ^ are almost entirely fixed by 
the field theory's chiral anomalies: 

b,,c , no rr,2 2n a ( L ,,b,,c,,d , no ..brr>2 i _/r>3 



& = CU^V + 2/3 a T^ - — ^- -C^ifif + 2/3 6/ i<T + 7 T 

-(-B) _ r< .fi ^ / n ,,C,,d I fl rri2 



Cb — Cabc/^ c : — I -zCbd/J, /J, + f3 b T 

e + P\2 J (5) 

I = ^CU^V + 2/3 a /i a T + 7 T 2 , 
S B) = -^WV + (3 a T . 
Here, C a & c is a symmetric tensor of anomaly coefficients: 

Cabc = ^ tr {T (a T fe T c) }, (6) 

where T a are the symmetry generators in the fermion representation. /3 a and 7 are numerical 
constants, which are not constrained by the anomaly. Note that a nonzero 7 is possible only 
for a parity-breaking theory. (3 a , on the other hand, is allowed in parity-conserving theories, 
as long as it is nonzero only for values of a corresponding to axial charges. Also, the group 
structure forbids a nonzero /3 a when a corresponds to a non-abelian charge. 



Our derivation generalizes the results of 



1] by allowing for gauged and/or non-abelian 



charges. We also elucidate the precise definition of the currents used in the derivation, 
leading to the correct value for Cqk~. Finally, we include the additional vorticity coefficients 
f3 a and 7, which were omitted in The results are consistent with our previous derivation 
for global non-abelian charges based on a null horizon dynamics in a gravitational dual 
description j^J. 

The appearance of the vorticity term in the anomalous hydrodynamic current has been 
first observed in the dual gravitational description of relativistic conformal hydrodynamics 
, |4j. Experimental manifestations of the vorticity term in heavy ion collisions have been 
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proposed in 

flfl 

. We note that hydrodynamics with SU(2) global symmetry charges has 
been considered in a holographic setup in p]]. 

The paper is organized as follows. In section [III we introduce external fields and carefully 
define the anomalous currents and their non-conservation equations. In section IIII[ we run 



through the derivation of eqs. ([3]). The derivation closely follows the one in jl|, with gauge- 
covariant derivatives introduced where necessary. It is essential for the argument that the 
currents and (non-)conservation equations defined in section [Til are gauge-covariant, and that 
C a bc assumes a symmetric form. Finally, in section IIVI we consider global currents whose 
conservation is slightly violated by anomalies with gauge fields. 

II. EXTERNAL FIELDS AND DEFINITION OF THE CURRENTS 

In order to better exploit the second law of thermodynamics, we will extend the theory 
by coupling the global hydrodynamic currents to external gauge fields. After this step, there 
is a gauge field A ^ coupled to each hydrodynamic current. The external gauge fields will 
enter along with the dynamical fields into the anomalous non- conservation equations. The 
external field strengths are of order ~ e, and are approximately homogeneous within each 
equilibration region; they may vary over the length scales of the hydrodynamic gradients. 
The covariant derivative will be understood to take the external fields into account. 

In a theory with anomalies, two subtleties arise with regard to the definition of the cur- 
rents and the precise form of their non-conservation equations. First, the anomaly diagram 
may be regularized in different ways, by adding different local counterterms to the action 
js]. These counterterms are given by local functionals of A®, and affect the definition of 
the current. Furthermore, the current may be defined not directly from the variation of the 
path integral, but again with the addition of some local functional of A a ^ [9| . Since there are 
no anomalies among the dynamically gauged currents, these issues only arise after we add 
the external gauge fields. What we would like to have is a set of currents that are gauge- 
covariant with respect to both the dynamical and the external gauge fields. The covariance 
can then be used to constrain the hydrodynamic constitutive relations. We will now see that 
there is essentially one such choice of currents, regardless of the regularization. 

Let the anomalous diagrams be regularized in any permissible way (i.e. keeping the dy- 
namical gauge symmetries intact). Consider a spacetime region Q, with boundary conditions 
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on dfl corresponding to our hydrodynamic state. Let us perform the path integral within Q 
with these boundary conditions, without integrating over the gauge fields: 

W[g^ u ,A°} = -zln / VipV$V$e iS ^ A »^^ , (7) 

where if) are the fermion fields responsible for the anomaly, and are any other fields present 
in the theory Our notation shows the dependence of W on the values of g^ u and within 
Q, and suppresses its dependence on the boundary conditions on dfl. The gauge fields 
include both the dynamical fields, whether or not they enter the fluid description, and the 
external fields. 

We now wish to derive the currents and the stress-energy density by varying W with 
respect to A^ and g^, respectively. However, in the hydrodynamic context we do not want 
the macroscopic currents to contain the D v F lLV contribution derived from the kinetic term 
of the dynamical gauge fields; also, we do not want the stress-energy density to contain the 
stress-energy of the macroscopic piece of the dynamical gauge fields. Therefore, we first 
subtract from W the kinetic term of the macroscopic gauge fields: 

W' = W + \J d'x^g- (F%) (FT) , (8) 

where the index a runs over the dynamical gauge fields, and the brackets denote averag- 
ing over the local equilibration region. This adjustment is gauge- invariant, so the gauge 
transformation properties of W are the same as for W . 
We now define the "consistent" anomalous currents as 

f a = SW'/SAl . (9) 

The divergence DfJ£ generates gauge transformations inside Q. The j£ are not gauge- 
covariant, but they can be made covariant by the addition of a local functional of A^. This 
is best explained in two steps. First, let G be the counterterm to the action which would be 
necessary to pass from our chosen regularization of the anomalous diagrams to the standard 
symmetric one (which would violate the dynamical gauge symmetries). Now, let us add 
SG/5A° to our currents. The resulting currents have the same transformation properties 
and satisfy the same non-conservation equations as the "consistent" currents under the 

n 

symmetric regularization. These are still not gauge- covariant. However, as shown in |9J, 
we can add to them an additional local functional which brings the currents into 
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a gauge-covariant form. We denote the new covariant currents by Jg. They satisfy the 
non-conservation equation: 



where C a ^ is given by (JBJ), with the indices a/?7 of the full symmetry algebra instead of 
abc. We stress that eq. ( fTUl) depends only on the matter content of the theory, and not on 
the choice of regularization. Also, in this context there is no conflict between the symmetric 
form of (1101) and the exact gauge invariance of the dynamically gauged currents. 

There are essentially no other choices of a gauge-covariant current: the only gauge- 
covariant functional of with the right dimension that can be added to J£ is D U F£ U . 
For the hydrodynamic currents, such a term was intentionally excluded. In any case, the 
divergence of this term vanishes identically, so it would not change eq. ( jTUl) . We also note 
that if one is interested in the constitutive relations for some other choice of currents, one 
should simply take the results (CQ) and (T5]) for J% and add to them the appropriate gauge- 
field functional. In principle, we should have another contribution to D^J^ from anomalies 
with gravity. However, these are proportional to the square of the Riemann tensor, which is 
0(e 4 ). The gravitational contribution is therefore negligible at the relevant hydrodynamic 
order. 

We now turn to the stress-energy density. So far, we have not performed the path integral 
over the dynamical gauge fields. Let us now integrate over all the gauge fields except for 
the macroscopic piece of the fields that participates in the hydrodynamics. This yields a 
new functional W'[g^ u , (Afy]. The stress-energy density is then defined as T^ u = 5W'/5g flu . 
This is a gauge-invariant quantity, since variations bg^ of the metric commute with gauge 
transformations 5A a , and the anomaly SW'/5A a is a metric-independent functional of gauge 
fields. From the transformation law under diffeomorphisms inside Q, we obtain the conser- 
vation law: 



n tm _ Z.CJ a e ^ vpa F 13 F 



(10) 




(11) 



where the covariant derivative on the RHS is defined in terms of \A°£j. Eq. (fTTj) can be 
written in a manifestly gauge-covariant way as: 




(12) 
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To see that the two equations are equivalent, we note that their difference is a functional of 
gauge fields alone (where the gauge fields that were already integrated out are understood 
as operators). Also, the difference must be gauge- invariant, since the LHS of both equations 
and the RHS of the second one are manifestly gauge- invariant. It must also be of mass 
dimension five. Finally, recalling the form of D v j v a and J v a — we see that the lowest 
power with which gauge fields appear can in the difference is three. There is no functional 
satisfying these conditions, so the difference between the two equations must vanish. 

It remains to complete the path integration, and to write the average of eq. (fTOj) for the 
hydrodynamic currents: 



On the LHS of this equation, the covariant derivative commutes with the averaging 
brackets, due to the weakness condition on the microscopic component of the dynamical 
gauge fields which participate in the hydrodynamics. By assumption (which will be relaxed 
in section HV)) . the RHS vanishes in the absence of external fields. In other words, it receives 
no contribution from products e M " pa F^F^. of two dynamical field strengths. Now, the 
external gauge fields do not have a microscopic component. Therefore, we can pull them 
out of the averaging brackets, giving: 



We've reduced the summed-over indices into the subspace of hydrodynamic currents, since 
only the gauge fields associated with such currents have a non-vanishing macroscopic com- 
ponent. 

III. CONSTITUTIVE RELATIONS 

The zeroth-order hydrodynamic constitutive relations are given by (TTJ. We will now use 
the non-conservation equations ffl2]) . f JT4|) and the second law of thermodynamics to derive 
the possible form of the first-order constitutive relations. From now on we will drop the 
averaging brackets, understanding that averaged quantities are always implied. Using the 




(13) 




(14) 
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definitions (jlj), we rewrite the relevant components of eqs. (fT2|) . (|T4j) as: 

(15) 

We fix the velocity w M as the unit timelike eigenvector of T" the energy density e as the 
eigenvalue of / yf—g corresponding to w' 4 , and the charge density as n a = u^J^ / y/—g. 
Using gauge covariance to constrain the possible contributions, the most general first-order 
corrections to the constitutive relations read: 



= - + V^h (£u»D v u v + xP^duP + ~a a P^D v ^ + a<&E 



(16) 



The various coefficient functions are constrained by the non-negativity of the entropy pro- 
duction rate. Using eqs. (115j) . the ideal fluid equations and the first law of thermodynamics, 
the entropy production rate can be written as: 



(17) 



The different terms in (1171) divide into those with a factor of e pvpu (via or B£) and 
those without. For the terms without & vpa , a standard exercise shows that we must have 
X = ( = } Xa = &a = &i E) = 0, cr^f ) = a ab , rj > 0, C > and a positive semi-definite cr (a6 ). 



a. 



The only way 



Let us now turn to the e^^-terms, closely following the derivation of 
for these terms to be non-negative is to vanish identically. The following consequences of 
the ideal fluid equations are useful: 

= 2a^ = —^(d^p - n a E«) , 

(18) 

P>iiP>a = -^E a „ + a,B p = -2^E a , - —B*(d& - n b E b ^) . 

We choose (p,/j,i/T) as the independent thermodynamic parameters on which the coeffi- 
cients (£ a , £„f \ £, fa 8 ^) may depend. We use the indices i, j, k to enumerate the hydrodynamic 
charges which are not screened, and can therefore assume nonzero values. For the gradients 
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of I and H B) in f JT7|) . we use the chain rule for the covariant derivative: 



where a is any thermodynamic function of state, with arbitrary charge indices. 
The condition on the e M ^ po -terms in f TTT|) can now be written as: 

T e + p y 

(20) 




(B) 



_i_ I 1 t(B) 1 r i , n «4 

At a point, the second-order terms on the right of each set of parentheses are independent 
Therefore, all the coefficients must vanish separately: 

oi 2f d£ ( a B) H B) 



dp e + p ' dp e + p 



(21) 



& JL = df (22) 



| + M = 2g B »; I^'f + ^ = V (23) 

T e + p T e + p T 

| = 2^ } ; ^ = ~<W (24) 

where the index q runs over the screened charges, for which n q = 0. To solve these equations, 
we use the thermodynamic identities: 

fdT\ T f dT \ rtT 2 



\d P ), i/T e + p \dOH/T)) p e + p' (25) 

Eq. fl2TI) then implies that £/T 2 and £a /T are functions of fii/T only. Using eq. (|22|) to 
eliminate & and £jjf , eq. f )23|) then gives: 

d(j/T*) _2 m d(ij B) /T) _ 1 

9(^/T) ~ T ' d^/T) ~ T ' 1 J 

For eq. (12"oT) to be integrable, we must have the symmetry conditions = C(y)& and 
Cijfc = C(ijk)- As we saw in section [TTJ, this is guaranteed by our construction of gauge- 
covariant currents. Integrating and taking into account eq. ( 1241 . we get the result (jSJ). The 
numbers (3 a and 7 arise as integration constants. Taking away the external fields, we obtain 
the constitutive relations of the physical fluid. 
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IV. NON-CONSERVATION DUE TO DYNAMICAL GAUGE FIELDS 



In the above discussion, we assumed that in the absence of external fields, the hydro- 
dynamic currents are exactly conserved. This assumption can be relaxed. To be part of 
the fluid description, a global charge only needs to be conserved approximately, so that its 
production rate is negligible within each local equilibration region. However, on the scale of 
hydrodynamic gradients, the conservation of the charge may be violated, as it was in sections 
HTllHII due to anomalies with external gauge fields. The same kind of "soft" non-conservation 
may take place in the physical fluid with no external fields, due to anomalies with the dy- 
namical fields. Taking into account both kinds of anomalies, the averaged non- conservation 
equation (TH1) takes the form: 

? (27) 
$ a = -C a B^ vpa (F^Fl) . , 

a g ap7 \ fxv pa/ mlcro ' 

where $ a is the charge production rate due to anomalies with the microscopic component of 
the dynamical gauge fields. Due to the group structure, <3> a can be nonzero only for global 
U(l) charges. Like every quantity in the hydrodynamic regime, $ a will be a local functional 
of the thermodynamic parameters and the macroscopic gauge fields. 

The $ a -term provides a correction to the hydrodynamic charge production rate in ( 1T5|) : 

D^Z = C abc E b ^ + $ a . (28) 

This propagates into the entropy production rate (fT7|) . which becomes: 

V = - (-T«^LV„ + J (1W (E^ - TDj±) - C ahc[ i a E\B^ - ^ 

In powers of the small parameter e, we are interested in $ a up to second order. A zeroth- 
order contribution to $ a must be suppressed by some small dimensionless factor, in order 
to maintain the approximate conservation of charge within each local equilibration region. 
Such a term does not mix with any other term in fT25]) . Therefore, its contribution to the 
entropy production must be always non-negative, i.e. /i a $^°^ a < 0. A first-order contribution 
$o also cannot mix with any other term in (|29|) . Such a contribution is then ruled out by 
the second law of thermodynamics, since it must be proportional to -D^m m , with no way to 
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(29) 



constrain its sign. Finally, a second-order term $a is not ruled out by the second law, and 
it will mix with the other terms in ( 129|) . In fact, such an arbitrary second-order contribution 
would affect the entire set of constraints on the transport coefficients, and invalidate our 

(2) 

conclusions. Thus, the only situation in which we can make a clear statement is when $a 
can be neglected, and only must be taken into account. In fact, this is a reasonable 
assumption: this is what we get if the same suppression factor present in also affects 

(2) 

$a , in addition to the suppression from small derivatives. 

In conclusion, we assume that $ a is given by some thermodynamic function of state, 
without derivatives or gauge field factors. This function must be suppressed by some small 
number, and satisfy // a $ a < 0. Our results for the transport coefficients remain unchanged, 
and the entropy production rate reads: 

V = ^ (2177^ + 3((DX) 2 + * ab (E a » - TDj±) (if - TD^)) 

1 (30) 
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